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Abstract. - We present and analyze a topologically induced transition from ordered, synchronized to dis- 
ordered dynamics in directed networks of oscillators. The analysis reveals where in the space of networks 
this transition occurs and its underlying mechanisms. If disordered, the dynamics of the units is precisely 
determined by the topology of the network and thus characteristic for it. We develop a method to predict 
the disordered dynamics from topology. The results suggest a new route towards understanding how the 
precise dynamics of the units of a directed network may encode information about its topology. 



Networks of interacting units prevail in a variety of systems, ranging from gene regulatory net- 
works and neural networks to food webs and the world wide web [1,2]. A fundamental question 
is: What kind of dynamics can we expect given a network of prescribed connection topology [3]? 
Even in networks of known dynamical units, known type of interactions between them and known 
topological details, it is hard to infer which kind of typical collective dynamics the network will dis- 
play (cf. Refs. [3-6]). If parts of the network exhibit residual symmetries, such as permutation or 
translation invariance, some general properties of the dynamics can be deduced [6], If, however, no 
symmetries remain, it is still an open question how topological factors can control network dynam- 
ics. See, e.g., [7-9] for some interesting recent approaches for phase oscillator networks of specific 
connectivities. 

In this Letter, we study directed networks of phase oscillators that exhibit a mechanism to synchro- 
nize and reveal general principles about how topology controls dynamics: Specifically, in networks 
with an invariant state of in-phase synchrony we analyze how the topology can control the units' dy- 
namics in a neighborhood of synchrony. We find that such networks, depending on their coarse-scale 
topological properties, belong to one of two classes exhibiting very different long term dynamics: Net- 
works of class I show in-phase synchrony in which each unit displays identical dynamics, independent 
of the unit's topological identity [10], i.e. independent of where in the network it is located. Networks 
of class II, instead of synchrony, show disordered dynamics. Here, together with the initial condition, 
the fine-scale topology precisely controls the dynamics of each unit. The dynamics therefore strongly 
depends on where the unit is located in the network - its topological identity. We develop a method to 
predict the disordered dynamics from the network's topology. Due to their topological origin, both the 
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Fig. 1 - Synchronization-disorder transition in directed network dynamics. (a),(b) The long time dynamics of two 
different random network realizations of N = 100 units with the same connection probability p = 0.05 started 
from the same random initial condition. The phase differences A<f>i := 4>i — ((f) j) . with respect to the average 
phase (4>j) j of all units are shown versus the units i, plotted relative to their possible range, 2n. 

separation of the ensemble of networks into two unique classes and the specific disordered dynamics 
realized by a network appear to be general phenomena and not restricted to the system studied here. 

Let us elaborate these findings. Consider a network of N phase oscillators i that interact via 
directed connections. The network topology is arbitrary and determined by the sets In(i) of those 
units j that have input connections to i, denoted j — ► i. We analyze a simple, paradigmatic model of 
interacting periodic oscillators, the Kuramoto model [11-13] defined by 

^-<pi(t)=Ui+ ^ Jij sin(0j - <f>i) (1) 

jeln(i) 

where the phase variable 4>i(t) e [0, 2tt) (with periodic boundary conditions) determines the state of 
unit i at time t, uji is its frequency, and Jij > is the strength of coupling from j to i with Jij = if 
there is no connection j — > i. In order to stress the topological effects, we neglect inhomogeneities in 
the dynamical parameters: we consider identical units = lo and homogeneous total input coupling 
strengths such that J^j Jij — J ( m a H illustrating examples we choose Jjj = J j fcj if unit i receives hi 
input connections from other units j). Without loss of generality, we take J = 1 in the following. We 
consider initial states in a neighborhood of the in-phase synchronous solution to reveal those features 
that apply for other oscillator networks as well. 

Observing the dynamics Q on different topologies, it was intriguing to find that seemingly sim- 
ilar networks (such as realizations of networks with identical degree distribution) yet displayed very 
different dynamics. Consider for instance the long-term dynamics of random networks in which ev- 
ery connection j — s> z is present with probability p. Several such networks show in-phase synchrony 
(cf. Fig. ^). Thus, the final states of the units display no information about the network topology. 
The units' topological identity is hidden. Other networks with identical statistical properties display 
disordered periodic dynamics, even when initialized in the same state (Fig[2)). In such disordered 
states almost every unit displays a different phase. It turns out (see below) that the network topology 
precisely controls the dynamics of these units. The units' dynamics thus display their topological 
identity! 

This phenomenon raises a number of questions. In which networks and how does the disordered 
state emerge? What determines the individual units' dynamics if the network displays disorder? 

To answer these questions, we studied the dynamics of small networks that exhibit qualitatively the 
same disordered state (see, e.g. Fig . |2^) . We tried to find a systematic dependence of the units' states 
on the network topology. As a first step, we ordered the units of the network (Fig. 2a) such that units 
with similar phases are displayed at proximate positions (Fig. 2b). This ordering reveals a division of 
the network in terms of its strongly connected components (SCCs) [14]. Whereas there are some units 
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Fig. 2 - (color online) Dynamics of small networks (N = 11) started from the same random initial condition. 
The upper column of each panel displays the directed networks with units labeled i £ {1, . . . , N}. The lower 
displays the relative phase differences A(j>i/(2n) versus i. (a) A network with homogeneous in-degree ki = 2 
for all units i exhibits an irregular asymptotic state, (b) Same network as in (a) with units with similar phases 
grouped, (c) A network with only one more directed edge (8 — > 6) [red; red dashed edge in (a)] compared to that 
of (b) induces a completely ordered state with identical dynamics of all units. 



with identical phases, several phases appear uniquely. In a network with just one more edge (Fig. 2c), 
the collective dynamics is completely synchronized such that the topological identity of all units is 
hidden (cf. also Fig. la). The above finding that ordering of the phases seems to reveal information 
about the coarse scale network topology led us to hypothesize that the partition of the network into 
SCCs is important to understand its dynamics. 

To test this hypothesis, we first analyze the dynamics of networks of arbitrary connectivity in a 
neighborhood of in-phase synchrony (</>,; (t) — <pQ (t) for all units i and all times t). Disconnected parts 
of a network can be treated independently, such that we focus on connected networks here. Sufficiently 
small perturbations Si(t) := 4>i(t) — <fro(t) to the synchronous state satisfy 

Sj = ^ J tj sm(Sj - Si) (2) 

for all i, which in first order approximation reads Si = — J Si + X^ein(i) Jij&ji or ^ = MS in matrix 
form, where 

( -J if j = i 

My = <^ j tj ifj e Hi) (3) 
[ ifj<£ {i} U In(i) 

are the matrix elements of M and S — (Si , . . . , Sn) t is the vector of the individual units' perturbations 
Si, This results in the first order period-T map (T = 2ir/ui) given by 

S(T) = AS(0) (4) 

where the matrix elements of A — e MT satisfy Ay > 0, reflecting the attractive couplings Jy > 0, 
and J2j Aij = 1 due to time translation invariance of the periodic orbit. 

For networks of arbitrary connectivities, this implies, via the Gersgorin disk theorem [15], that 
all eigenvalues A; of A satisfy |Aj| < 1. A sufficiently small perturbation to the synchronous state 
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Fig. 3 - (color online) Decomposition of networks shown in Fig. 2 in terms of their strongly connected components 
(SCCs). On the left, the vertices of the networks are grouped to SCCs s 6 {1, 2, 3, 4} (large italic numbers). On 
the right, the level structure of these components is shown, (a) Three-level network with two source, cf. Fig. 2a,b. 
(b) One additional link makes it a four-level network with one source, cf. Fig. 2c. 

cannot grow (in maximum norm), cf. [5], such that synchrony is at least marginally stable. Moreover, 
independent of the network connectivity there is one eigenvalue Ai = 1 with an eigenvector vi = 
(1, 1, . . . , 1) corresponding to the uniform phase shift. 

If the network is strongly connected [14] the Perron-Frobenius theorem [15] guarantees that the 
largest eigenvalue Ai = 1 is unique and all other eigenvalues satisfy |A,| < 1 fori £ {2, . . . , N}. This 
implies that the synchronized state is asymptotically stable and thus locally attracting. In networks 
of irregular topology we even often find that the system converges towards it from arbitrary initial 
conditions. 

If the network is not strongly connected it consists of two or more strongly connected components 
(SCCs) and the analysis of the asymptotic dynamics is more involved. For better accessibility of the 
main points of this Letter we describe the details of this analysis in the Appendix. Briefly, for a given 
network, we first determine the SCCs and the uni-directional connections among them. Second, we 
determine the level structure of this super-network of SCCs (cf. Fig. [3}. Based on this composition 
analysis we have revealed a number of distinctive features of the dynamics on directed networks: 

The ensemble of networks divides into two classes with qualitatively different long term dynamics 
(initialized sufficiently close to the in-phase solution). All networks that have M = 1 source SCC 
(which does not receive any input from other SCCs) belong to class I: this source SCC is guaranteed 
to synchronize because it itself is a strongly connected network without further input. Since each unit 
i performs a local weighted averaging of phases determined by the weights Ay in all units outside 
the source component asymptotically converge towards the (common) phase of the units within the 
(only) source SCC. This result also follows explicitely from the analysis given in the appendix in the 
special case of only one source component in level I = 1 (and no source components in levels I > 1). 
It implies that for all networks with one source SCC the local asymptotic dynamics is also in-phase 
synchrony. In contrast, networks having M > 2 source SCCs (class II) typically show disordered 
dynamics. These M source SCCs can synchronize independently of each other, creating M — 1 
independent phase differences which result in an {M — 1) -dimensional continuous family of periodic 



Marc Timme: Does dynamics reflect topology in directed networks? 



5 



XXX 

■ • X 

X 

. . . 

XXX 



71 68429113510 
unit i 

Fig. 4 - Prediction of the dynamics in a disordered state (cf. Fig. 2a,b) based on the composition analysis. The 
relative phase difference A<£j/(27r) is shown for the grouped units i. The linear prediction (x) of the actual 
phases (•) (based on one intial state) well distinguishes the ordered from the disordered state (which would 
be a constant at zero) and even is a good indicator of the quantitative dynamics of the units. The asymptotic 
phase dynamics started from a different initial state (gray •) illustrates that in this example other initial states 
yield a pattern that is distinguished from the former pattern only by a real mulitplicative factor (in first order 
approximation). 

orbits, that include the synchronous state as only one specific orbit. All these orbits are marginally 
stable, in particular the synchronous state has a basin of attraction of measure zero, such that the 
dynamics is almost surely disordered. For the examples above, we find that the dynamics shown in 
Fig. la originates from a class I network whereas that of Fig. lb originates from a class II network. 

The composition analysis also reveals how the details of the topology of the network precisely 
control its dynamics in the disordered state: The topological identity of each unit, particularly the 
fine scale topology of that SCC it is part of, determines the unit's dynamics. In fact, we can predict 
the disordered dynamics on a fine scale: Given the initial state 0(0), we uniquely determine the 
approximate phases of all units recursively level by level, and hence predict the complete collective 
dynamics of the network from the topological identity of their units (see Appendix). Figure0]illustrates 
such a prediction. It resembles well the actual dynamics of the units. 

Reversely, partial information about the topology of the network may be obtained from knowing 
the disordered dynamics of its units: Iterating Eqn. (|5}, we obtain explicit linear restrictions of the 
space of all networks from the disordered dynamics by imposing its invariance. So only a lower- 
dimensional subset of networks is consistent with the phase pattern. 

What is the mechanism underlying the transition to topology-induced disorder? The following 
description is general; nevertheless it is instructive to imagine, as an illustrating example, a network 
composed of two source SCCs and one sink SCC which receives input from the other two. A strongly 
connected network, and thus each source SCC, synchronizes completely. However, different of these 
source SCCs typically converge towards different phases, that depend on the initial state. If now 
different units in a downstream SCC are pulled towards different phases, and there is a complicated 
pattern of connections between them within this SCC, the dynamics of all its units will typically be 
distinct. In particular, the units' dynamics depend on the phases of the units in connected upstream 
SCCs, i.e. indirectly on the initial state of the network and on the specific topology of the SCC 
considered. 

All these phenomena appear to be general and are not restricted to the model system Q considered 
here. This is due to the topological origin of the phenomena: First, the transition line between networks 
of classes I and II is identical for various kinds of oscillator networks having an invariant in-phase 
solution. Second, given an initial state sufficiently close to synchrony, the disorder in the long term 
dynamics is characteristic for the topology of a network: The linear analysis (see Appendix) holds as 
well for all disordered dynamics that are topologically equivalent within the class of periodic single- 
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variable oscillator networks. We checked this explicitely for networks of (i) Kuramoto oscillators 
Q with coupling functions different from the sine function and (ii) spiking neural oscillators where 
interactions are delayed and mediated by pulses that occur only at discrete instances of time [5, 17]. 
Although we have no proof of how general these results are beyond single-variable oscillators, we also 
observed that even (iii) networks of diffusively coupled chaotic Rossler systems [18] behave similarly. 
On the same network topology, the dynamics of these three kinds of distinct systems show closely 
related patterns of phase disorder. 

Commonly, transitions from synchrony to disordered dynamics have been devoted to hetero- 
geneities of, e.g., dynamical parameters or the degree distribution, cf. [19-21]. However, the pre- 
cise impact of topology onto the dynamics of directed networks revealed here was so far not noticed. 
Even recent studies, considering the exact dynamics of networks of given, specific topologies (see, 
e.g. [5,22] and refs. therein) have not taken notice of this impact. The main reason for this may be that 
all example networks chosen to look at explicitely again were standard cases such as highly connected 
random networks or lattices. 

Real-world oscillator networks, occurring across disciplines in physics, biology and technology 
[3, 13], however, have a far more complicated topology, and, as demonstrated in this Letter, may 
thus strongly deviate in their dynamics. In a study [23] related to ours, the notion of long range 
action has been introduced showing that in certain directed networks of iterated maps the dynamics of 
boundary units can control the dynamics of the entire network. Our results suggest that local and global 
topological features, such as the SCC super-network and the detailed topology of particular SCCs, 
may act together to precisely control the dynamics of individual units in complex directed networks. 
The concepts developed here may thus also help to uncover information about the topologies of such 
networks from their dynamics. 

I thank M. Denker, M. Diesmann, T. Geisel, C. Kirst, A. Levina, R.-M. Memmesheimer, F. Wolf, 
R Ashwin, L. Bunimovich, J. Borresen, S. GroBkinsky, B. Kriener and S. Strogatz for stimulating 
discussions and comments on the manuscript. I acknowledge partial support by the Federal Ministry 
of Education and Research (BMBF), Germany, under grant number 01GQ0430. 

Appendix: Level Structure of the SCC Super-Network: The level structure of the SCC super- 
network is constructed in three steps. First, we determine the SCCs of the network using a standard 
method [24], the computational complexity of which is O(N). Second, connections between them are 
straightforwardly derived from the underlying connections between units comprising these SCCs. A 
connection from one SCC to another, s — > s' is present if there are i E s and j £ s'with a connection, 
i — > j between them. Third, we find the longest undirected path from any source SCC (without 
incoming connections) to any sink SCC (without outgoing connections). The length of such a path is 
found by counting a connection followed along its direction as "+1" and against its direction as 
All units i in a source SCC of the longest path is given the level number £(i) = 1. The levels of all 
other SCCs are determined recursively according to the above counting rule. The computational costs 
of finding the inter-SCC connections and the level structure strongly depend on the network under 
consideration. 

Dynamics From Network Topology: Given the level structure, the linearized dynamics of every unit 
is determined for all units in every given level, starting with level £ — 1 and proceeding through sub- 
sequent levels recursively. Let 4> = (rf*^ > • ■ • > </>^J = , ■ ■ ■ , fa N ) denote the asymptotic phases 
of all units fa in terms of the collection of phases <f)^' of the units at a given level £ £ {1, . . . , L}. 
For all units i with £(i) — 1, their final states are 4>f^ — c s , where c s depends on the initial state 
4> s (0) restricted to the SCC s. It equals the first component of the vector c s = V~ 1 4> s (0), where 
V = (vjj , . . . , v^) is a matrix of the R eigenvectors v ir of A localized on the SCC s and v il is the 
eigenvector corresponding to the eigenvalue =1. This yields the vector (fr- 1 ' of asymptotic phases 
in all units in level t = 1. 
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The phases <jr ' of units in the other levels I > 2 are determined iteratively given the phases cjr ~ ' 
of units in level I — 1. If some of the <f>^ are part of a source SCC in level £, these phases 4>$ UIce 
are determined analogous to those in level 1=1, The corresponding sub-matrices A^_i and Ae£ of 
the matrix A in (0) needed to determine the remaining phases SO urce describe the interactions with 
units of the previous level £ — 1, and within the SCCs of the current level £, respectively. Note that by 
definition of the level structure, there are no interactions from level I to level I — 1. Thus the equation 
encoding this uni-directional dependence, ^source = 4e r e-i0^~ 1) + At,t<$L Source- yields 

Source = (1 - A^A^.^-V (5) 

such that, together with the Source fr° m above, all phases <p^ of units in level £ are determined. 
Iterating this for all levels £ £ {2, . . . , L) we obtain the linear prediction of the complete disordered 
asymptotic state <j>. This analysis only depends on the linearized effective couplings Aij that determine 
the SCC super-network and applies thus not restricted to the system Q. 
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